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Abstract 

We examine time ordering effects in strongly, suddenly perturbed two-state quantum systems 
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(kicked qubits) by comparing results with time ordering to results without time ordering. Simple 



^i analytic expressions are given for state occupation amplitudes and probabilities for singly and 

multiply kicked qubits. We investigate the limit of no time ordering, which can differ in different 
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I. INTRODUCTION 

There are two reasons to consider time ordering in kicked qubits. First, the behavior of a 
two-state quantum system interacting with a rapidly changing external field, i.e. a diabat- 
ically changing qubit, may be described analytically. While such solutions were examined 
some forty years ago in the context of Landau- Zener transitions in atomic and molecular 
reactions lit], relatively little attention has been paid to this problem in the context of more 
recent work using two-state systems [2| ranging from quantum computing |3|, |4| to quantum 
control of atomic and molecular reactions 5] to manipulation of matter waves |6{], where this 
class of analytic solutions may be useful. The second reason is that time ordering has been 
used recently p, la, |9| to formulate an understanding of time correlation in multi-particle 
systems (or, in the context of this paper, systems of interacting qubits). The central question 
here is how one particle (or qubit) is connected with other particles (or qubits) in the time 
domain. This problem has previously been formulated J8| using second order perturbation 
theory, where observable time correlations between different particles arise from time order- 
ing of weak, external interactions in atomic scattering [lOj . The kicked qubit gives us an 
opportunity to study the nature of time ordering for a simple, analytically tractable system 
in a strongly non-perturbative regime. 

Except for relatively simple e~ tEjt / h phases, where the Ej are simple eigenvalues of a time 
independent Hamiltonian, there are only two ways in which time enters the time evolution 
of a quantum system. The first is through the explicit time dependence of an external 
interaction V(t), and the second is through the constraint of time ordering imposed by the 
time dependent Schrodinger equation itself. This time ordering, discussed below, imposes a 

causal-like constraint that places operators such as H(t n ) H(i,2)H{ti) in order of increasing 

time. This confining condition interrelates the influence of the time parameters t n . . .t\. In 
second order perturbation theory it has been shown JlOJ] that the time ordering constraint has 
negligible effect if either V(t) or its variation with time is sufficiently small over the time of 
the experiment (perturbative or constant potential limit) or if the energy levels of the system 
before perturbation are all nearly the same (degeneracy limit). In either case, principal value 
contributions from energy fluctuations in short-lived intermediate states vanish (111]. 

In this paper we formulate the problem of time ordering in a non-perturbative two- 



state quantum system, i.e. a qubit. After describing the basic formalism in Section III A^ 



in Section III Bl we define the limit without time ordering and show that time ordering 
disappears in either the constant potential limit or in the limit of degeneracy of the two 
unperturbed states, where H(t') and H(t") commute. We discuss the relationship between 
time ordering and the adiabatic approximation. Analytic solutions for singly and multiply 
kicked qubits are presented in Section Hi CI with and without time ordering. In the case of 
a single kick, we show how time ordering affects the transition probability from one state to 
another. In the subsequently discussed case of a double kick, any transition is due entirely to 
time ordering effects. We discuss corrections for pulses of finite duration and in Section HTTI 
provide calculations illustrating our results. We present most calculations and some key 
formulas in both the Schrodinger and intermediate (or interaction) pictures, and discuss 
some differences between time ordering effects in the two pictures. 



II. THEORY 



A. Basic formulation 



Consider a two-state system, whose states are coupled by a time dependent external 
interaction, e.g. a qubit with "on" and "off" states. The time dependent Hamiltonian for 
this system may be expressed as 



H(t) --■■ H + V(t) 

-AE/2 
AE/2 



+ 



AE 



cr z + V(t)a x 



V(t) 
V(t) 



(1) 



where AE = E 2 — E\ is the energy difference of the eigenstates of Hq. Here a x and a z are 
the usual Pauli spin matrices. 

Two simplifying, but removable, assumptions have been made in the second line. First, 
we assume that all of the time dependence in the interaction operator V(t) is contained in 
a single real function of t, which is often justifiable on experimental grounds 12|, |l3|, Il4j . 
Secondly, in this paper we assume for convenience that the interaction does not contain 
a term proportional to H . Obviously, an interaction operator V having the form of a 
combination of a x and a y (i.e. a complex time dependent field) is equivalent to the above 



form after rotation of coordinates. As we shall discuss later, there are other choices of how 
to separate H into Hq + V, and these choices have consequences. 



For a qubit with "on" and "off" states 
evolve according to 
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the probability amplitudes then 



ih 



d 

It 



a i 

«2 



-AE/2 V(t) 
V(t) AE/2 



Oi(t) 

a 2 (t) 



(2) 



The solution to Eq. (J2J) may be written in terms of the time evolution matrix U(t) as 



CLi(t) 

a 2 (t) 



U(t) 



01 (0) 

02 (0) 



U n {t) U 12 {t) 
U 21 (t) U 22 (t) 



oi(0) 
o 2 (0) 



(3) 



where an experiment is begun at a time t — and completed at t — Tf, Since we assume 
the two-state system is closed, P\{t) + P 2 (t) = |ai(t)| 2 + |a2(t)| 2 = 1. 
The time evolution operator U(t) may be expressed here as 



u(t) = T e -H *(*')<**' = T e -££(-¥^+n*>*K 

-i/h) n rt 



(4) 



n=0 



/?! 



H(t n )dt n . 



H(t 2 )dt 2 / Hit^dh 



The only non-trivial time dependence in [/(£) arises from time dependent -£/"(£) and time 
ordering T. The Dyson time ordering operator T specifies that H(tj)H{tj) is properly 
ordered: 

TH^Hfa) = H(U)H{t ) + 0^ - t t ) [H(tj), H(U 

Time ordering imposes a connection between the effects of H(ti) and H(tj) and leads to 



observable, non-local, time ordering effects [12 , 



uQ 



when 



H{h),Hlti) ^°- 



i. Pulses 



In this paper we regard V(£) as having the form of a smoothly varying pulse, or sequence 
of pulses, each of duration r and peaked at T&. We define phase angles 



n I V(t')dt'/h 
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f3 = rAE/2h 
jt = tAE/2h. 



(5) 



The angle a is a measure of the strength of the interaction V(t) over the duration of a given 
pulse. In this paper we are mostly interested in the non-perturbative regime corresponding 
to a > 1, so that substantial changes in the state occupation probabilities, Pi and P2, may 
occur. The angle /3 is a measure of the influence of H during the interaction interval r. 
The angle jt is the phase accumulation of the propagation due to H over a time t. The 
diabatic (kicked) limit corresponds to (3 <C 1, the perturbative limit corresponds to a < 1, 
and the adiabatic (slow) limit generally corresponds to r — ► 00. 

B. Time ordering 



Since time ordering effects can be defined as the difference between a result with time 
ordering and the corresponding result in the limit of no time ordering, it is useful to specify 
carefully the limit without time ordering. Removing time ordering corresponds to replacing 
T — > 1 in Eq. (jlj. This corresponds to the zeroth order term in an eikonal-like, Magnus 
expansion in commutator terms pj|. In the limit of no time ordering, a multi-particle time 
evolution operator factorizes into a product of single-particle evolution operators 



1 . Limit of no time ordering 



Replacing T with 1 in Eq. (jljl. in the Schrodinger picture we have, 

-i/h) n 



(J(t) = Te-Uo Ht')<tt' ^y, 

-i/h) 



n=0 



III 



H(t')dt' 



n=0 



IV. 



H t+ / V(t')dt' 





n 00 
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n=0 



-i/h) n 



IV. 



e -iHt/K = (j^ t) 



(Ho + V) t 



(6) 



where 



Vt = [ V(t')dt' = f V{t')dt'a x = aho-z 
Jo Jo 



H = Hq + V, and [Ho, V] terms are non-zero. By expanding in powers of [H(t"), H(t')], it is 
straightforward to show that to leading order in V and H the time ordering effect is given 
by 

1 



U-U° ~ ? / dt" / dt' 

2h 2 Jo 



H(t"),H(t>] 



2fr 



H ,V / dt'(t-2t')f(t'), (7) 
J Jo 



where V{t') = V f(t'). This leading term disappears if the pulse centroid T k = t/2 and f(t') 
is symmetric about T k . Furthermore, U — U° vanishes identically in the special cases of 
V(t') = 0, V(t') = V, or AE = 0, as will be discussed below in Section lHB2l \l(\. 

In general there is no simple analytic form for the exact result U(t). For the result without 
time ordering, we have 



U°(t) 



-Hyta z —iaa x 



cos Z + ijt^ 



—la 



sing 



-ia-JRi cos ^_i 7t M|l 



(8) 



where £ = J a 2 + ijt) 2 . Here we have used the well known identity |l7( e l ^ a ' n = lcos0 + 
iff ■ h sin <j), following from (ff ■ h) n = 1 (or a ■ h) if n is even (or odd). 

Similarly, in the intermediate, or interaction, picture, Ui(t) = e lHot U(t), and one has 



Uj{t) = Te-tti*™* 

{-i/ny 



JIM*)** 
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n=0 



n\ 



-in oo 



Viffldt? 
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n=0 



-l/h) n r 



IV. 



Vit =U%t) 



(9) 



where V}(t') = e iAot ' l % V (t')e- i&ot ' l % and V it = / * Vi(t')dt'. For a Gaussian pulse of the form 
discussed in Section ITTTl 



Vit = ahe & [a x cos 2^T k + & y sin 2^T k \ 



(10) 



and 



U°(t) 



-2vyT k 



(11) 



cos f ae ^ J — isinlae" 

-isin(ae _/3 ) e 2%lTk cos iae~^ ) 

as long as the measurement time t is after the completion of the pulse, i.e. t — T k ^> r. 

It has been shown previously [lOj that to second order in perturbation theory, II i — Uf ~ 
Vi(t"), Vi(t') , somewhat similar to the commutator in the Schrodinger picture above. From 
this we immediately see that time ordering effects do not appear until second order in a 
perturbative expansion in a. Again, Uj — Uj — > in the special limits Vj — > 0, Vj — > Vi, or 
AE — > 0, to be discussed immediately below. However, we will also find in Section Till B II 
that tli — Uj vanishes in the diabatic limit of a single ideal kick, V(t) ~ S(t — T k ), whereas 
U — U° is non-zero. Thus, in principle, the definition of the limit of no time ordering depends 
on the picture (representation) used. In the discussion we shall relate this difference to the 



gauge choice of how one separates H into H + V(t). As shown in calculations presented 
below, this difference can be negligibly small under some conditions. 

In both pictures, time ordering effects are associated with the fluctuation of a time de- 
pendent interaction about its time averaged value. 

2. Relation to other limiting cases 

Now we compare the limit without time ordering with the degenerate, weakly varying 
potential, and adiabatic limits. The connection with the diabatic (kicked) limit appears in 
Section Hi C| where we discuss the analytic solution for U in the case of a short pulse. 

We noted above that for a general pulse, there is no analytic solution for U(t). However, 
in the limit when the unperturbed states become nearly degenerate, i.e. AE <C h/Tf, we 
obtain 



—iaa x 



sin rv 

(12) 



cos a —ism. a 
-i sin a cos a 



This result may be obtained [18[ either from the coupled state equations of Eq. (J2J, or from 
Eq. (J3J with AE —>■ 0. In this degenerate limit the mathematical complexity of the qubit 
simplifies significantly, as may be seen by comparing Eq. (J12JI with Eqs. f|15|) and (|17j). Most 
of the complex time connections have been removed. We call this degenerate qubit a dit. 
In such a dit if the phase angle a equals tt/2, then an "on" state is turned off and an "off' 
state is turned on with probability 1. Then the dit is further reduced in complexity to a 
trivial classical bit. 

We notice that Eq. (fT2*|) may also be obtained from Eq. (jHJ) by taking 7 — » 0. Thus, time 
ordering effects vanish in the degenerate limit. 

A second situation in which time ordering effects generically become small is the case of 
a constant or weakly varying potential, V — V <C h/Tf. Then the full evolution matrix 
approaches the non-time-ordered expression given by Eq. (JHJ or Eq. 0. Clearly the per- 
turbative limit a C 1 is a special case of this, but U — U° vanishes also for a large average 
external potential V, as long as fluctuations around V are small. Such a situation may 
sometimes be addressed more transparently by absorbing the average part of V into Hq. In 



the extreme case V = 0, Eq. (jHJ) reduces to 

U°(t) -► e 



iHot/h —iyta z 



e"^ 4 



(13) 



which of course agrees with the exact evolution matrix U(t). 

In summary, time ordering effects disappear either when (i) a < 1 or when (ii) /3 <C 1 
('jTf <C 1) in the intermediate (Schrodinger) picture. The physics becomes especially simple 
in the overlap of regimes (i) and (ii). There, one easily finds that in the Schrodinger picture, 

for example, 

1 + ijt —ia 

—ia 1 — ijt 

It can be shown that corrections and time ordering effects start at O (a^t 2 ) and O {a 2, yt) 
for a symmetric pulse centered at T^ = t/2 (see Eq. (JHU)). The situation is similar in the 
intermediate picture, except that time ordering effects vanish identically at leading order in 
a, and begin at O (a 2r yt) only. 



U(t) ~ U°(t) 



(14) 



H ,V 



0; however this 



Strictly speaking, time ordering effects also vanish if (iii) 
situation is of little practical interest due to the fact that no transition or population transfer 
is possible. 

In addition to the degenerate, perturbative, and diabatic regimes, a fourth limit exists in 
which analytic solutions for U(t) are generally available. In the adiabatic limit where the 
external interaction V(t) changes slowly in time [lflj . it is useful to define the instantaneous 
level splitting Q(t) = ^J(AE) 2 + 4V 2 (t) and the accumulated phase 9 = J 4 tt{t')dt'/2h. Then 

U(t) -> U\t) ^ (15) 

cos 9 cos 0_ + i sin 9 cos <ft + cos 9 sin </>_ — i sin 9 sin <; 
— cos 9 sin d)_ — i sin 9 sin <$> . cos 9 cos <b_ — i sin 9 cos < 



Here 0± = (<j)(t) ± 0(0)) /2 where <f)(t') 
V(t) = V(0), relevant for a pulse, 



tan : (2V(t')/AE). In the special case when 



U A (t) 



cos 9(t)+i^ sin 9(t) -2i^sin9(t) 
-2i^sin9(t) cos9(t) - i£^sm9(t) 



(16) 



We note that this solution is similar mathematically to the rotating wave approximation 
(RWA), which is widely used to describe atomic transitions using external fields tuned to 
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frequencies near the resonant transition frequency between two states |!9ll20Ll2ll|. Inserting 
Eq. (JT5|l into Eq. (J2J), one finds that the leading correction is small when hV(t')AE <C 
Q 3 (t f ). When the splitting AE of the unperturbed qubit is not small, i.e. AE > V, 
then this adiabatic validity condition reduces to the Landau-Zener criterion |22J, namely, 
hV(t') <C (AE) 2 . The leading correction to Eq. (J15J1 or Eq. (|16|) is then given by non- 
adiabatic transitions at the avoided level crossings, where V(t') and thus the level splitting 
Q(t') goes through a minimum. More generally, for a pulse having a smooth shape, such 
as the Gaussian pulses discussed in Section HTT1 the criterion %V(t')AE <C Q 3 (t') reduces to 
the union of a <C j3 2 and (3 <C a 2 . 

Remarkably, the adiabatic regime overlaps both with the degenerate limit (when (5, jTf — > 
at fixed a) and with the perturbative limit (when a — » at fixed /3, jTf). Thus, Eq. (I15|) 
reduces either to Eq. (|T2*|) or Eq. (fTSjl as AE — *• or V — > 0, respectively. In these overlap 
regions, time ordering effects are small. More generally, however, time ordering effects in the 
purely adiabatic regime (a ^> 1 and /3 > 1 or j3 ^> 1 and a > 1) are large. Time ordering 
effects are also large when a and (3 are both of order unity, where no simple analytic solutions 
exist for the full evolution matrix U(t). 



C. Kicked qubits 



Now we consider time ordering in kicked qubits 



, i.e. the diabatic limit where V(t) 



5(t — Tfc). First, we present analytic expressions [23|, |2J] for a kicked qubit, i.e. a two 
state system subject to an external interaction, V(t), that changes rapidly with respect to 
Tae = tt/7 = 2irh/AE, the period of oscillation of the free system. The corrections needed 
for finite pulses will be briefly analyzed. We also discuss the extension to multiple kicks, 
using the double kick as an example. Finally, we consider the influence of time ordering in 
these kicked systems. 

1. Single kick 

Here we consider a two-state system where the interaction V(t) may be expressed as a 
sudden "kick" at t — Tf., namely V(t) = ahS(t — T&). For such a kick the integration over 



time is trivial and the time evolution matrix in Eq. (@J) becomes 



u K (t) 



ei ^ TM n e -ij T *:iv{t')M« x /n ei ±M Tk(Tz/n 



Jj(t-T k ) 





-if(t-T h ) 



e' 7 * cos a 



-ze" i7(t_2Tfc) 



sin a 



cos a —i sin a 
—i sin a cos a 
j e n(*-2Tfc) s in a 



e * 7 * cos a 



3 *7 r fe 



-i7Tfc 



(17) 



for t > Tfc. The second line follows from the identity given below Eq. (jHJ) above. As explained 
below this solution is valid when (5 -C 1 so that there is little effect from ff during the short 
time when V is active. 

From Eqs. (J3J) and (J17J1 we have for a kicked qubit initially found in state 1: 



Pi(t) 



Mi)| 2 
Mt)| 2 



l^nWP 



cos a 



\um\ 2 - - 



sin a 



(18) 



The corresponding probabilities for a kicked qubit without time ordering are discussed below 
in Section 111 (Ml 



2. Finite pulse corrections for a single pulse 

When the pulse width is finite the corrections to Eq. (JT7|) are 0({3) and corrections to 
Eq. (J18J1 are 0(f3 2 ). These corrections result from the commutator of the free Hamiltonian 
Ho with the interaction V during the time r when the pulse is active. This is related to the 
series expansion that arises in the split operator method 25] . For example, in the case of a 
rectangular pulse of width r, the exact time evolution is given by 



/Vrectangular / j_\ 



^■yt—i/3 



j r i'y(t-2Tk) n . sma' 



(cosa'H-i/3^ 



(19) 



where a' = \fa 2 + f3 2 . To leading order in j3, i.e. in the width of the pulse, the error in the 
kicked approximation is given by 



SU{t) = Tjrect^gul^^ _ jjK ^ = -^ / 



/sin a 



— cos a 



a 



e i7i 
-e~ i7< 



(20) 
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For a narrow pulse having a generic symmetric shape, the leading correction to the kicked 
approximation will still have the form 



5U(t) = U(t) - U K (t) = i/3g(a) 



-e"* 7 * 



(21) 



where g(a) is now a function that depends on the shape of the pulse. By comparing Eqs. (J3J) 
and (J 17)) at leading order in AE, after some algebra one obtains 

2 



g( a ) = - fdt COS 2 (f V{t')dt'/h) - ros'in/2) 



(22) 



Expanding U(t) of Eq. (J3J) and U K (t) of Eq. (|T7|) simultaneously in AE and V, or equivalently 
in (5 and a, we find 



<JC/(t) = U(t)-U K (t) 



+ 



iA£ 



i(A£) ; 



dt 



t x 2 



v(t')dt'/h 



Jv(t')(t'-T k ) 2 dt'/h 3 



e i7t 
-e"* 7 ' 

q e vy(t-2Tk) 

e -h(t-2T k ) q 



(23) 



so the two leading correction terms scale as /3a 2 and j3 2 a. 

3. Multiple kicks 

A series of either identical or non-identical pulses can easily be handled by multiplication 
of several matrices of the form of Eq. (J17J1 . For example, one may consider a sequence of 
two kicks of opposite sign at times t — T\ and t = T 2 , namely, Vkick antikick(^) = ceK6(t — 
Xi) — aHS(t — T 2 ). Following the procedure given in Eq. (|17j) one obtains the time evolution 
matrix for t > T 2 , 



rkick antikick/ 



i^l{t-T2)cr z iaa x i^{T 2 -T{)a z -iaa x i^T-ia z 



t thick antiKicK/j\ 

e^ (cos 'jTs + i sin 7T S cos 2a) e n(*-2T) g j n y^ g j n 2 a 



(24) 



-i 7 (i-2T ) , 



sin *yT s sin 2a e l< > (cos 7TI, — i sin ^T s cos 2a) 



where ( = 1 (t-T s ),T = (T\ + T 2 )/2, and T S = T 2 -T 1 . As 7T S -> 0, ?7 kick antikick (£) reduces 
to Eq. (H3|). 
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For a double kick with V = 0, we have from Eqs. (J3J) and 

P x (t) = \a 1 (t)\ 2 = \U k { ck antikick (t)| 2 = cos 2 yT s + sin 2 yT s cos 2 2a 

P 2 (i) = |a 2 (i)| 2 = \U^ ck antikick (t)| 2 = sin 2 yT s sin 2 2a . (25) 

The single-kick result of Eq. (|17|) remains valid for two or more kicks of combined strength 
ax + . . . + a n — a if the total phase associated with the inter- kick free evolution, 7(T n — Ti), 
is small. A mathematical analysis for multiple kicks separated by arbitrary time intervals is 
straightforward, but not included here. In the case of a periodic series of pulses with period 
T, the time evolution may be obtained by diagonalizing the matrix of Eq. (|T7jl and finding 
the Floquet eigenstates and eigenphases. The two Floquet eigenphases are then given by 
e ±tx , where x — cos -1 [cos a cos 7T]. 

4- Time ordering for single and multiple kicks or pulses 

We first consider the case of a single kick or pulse. In the Schrodinger picture, time 
ordering effects are present even for a single ideal kick, specifically the time ordering between 
the interaction and the free evolution preceding and following the kick. Thus, in the absence 
of time ordering, the time evolution U°(t) is given by Eq. (JHJ), which differs from the exact 
expression U K (t) of Eq. (|17|) when a and yt are both non-zero. The time ordering effect 
t) (t) — U°(t) vanishes in either the degenerate limit yt —> or in the perturbative limit 
a — > 0. For small a and yt and assuming T^ = t/2, the time ordering effect at leading order 
takes the form of a sum of 0(a(7^) 2 ) and O (a 2 yt) terms. For a Gaussian-shaped pulse, 
the transition probability in the Schrodinger picture without time ordering is given by the 
second line of Eq. (|2~9|l . 

In the intermediate picture, time evolution without time ordering for an ideal kick is 
obtained by substituting (3 = into Eq. (|11|) . and agrees perfectly with the exact expression 
of Eq. (|17jl . when the latter is transformed into the intermediate picture. Thus, time ordering 
effects disappear for a single ideal kick in the intermediate picture, in contrast with the 
Schrodinger case. This is easily understood by considering that in the intermediate picture, 
time ordering is only between interactions at different times, V}(t') and V}(£"), not between 
the interaction V(t') and the free Hamiltonian H (t"), as in the Schrodinger case. For a 
single ideal kick, all the interaction occurs at one instant, and no ordering is needed. Of 
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course, for a finite-width pulse, i.e. j3 ^ 0, time-ordering effects do begin to appear even in 
the intermediate picture. To leading order, Uf-{t) — Uf(t) = 0(a 2 j3). We note that the time 
ordering effect in the intermediate picture is independent of the measurement time t, though 
it does depend on the pulse width r through the (3 parameter. For a Gaussian-shaped pulse, 
the transition probability in the intermediate picture without time ordering is given by the 
third line of Eq. (|25|). 

We are now ready to examine the time ordering effect for a multi-pulse sequence, fo- 
cusing on the pulse-antipulse scenario of Section III C 31 In the limit of no time order- 
ing one has P%{t) = in the Schrodinger picture as seen from Eqs. (JHJ) and 1)13)1 with 
V = 0. In the intermediate picture, however, the transition probability is non-zero even 
without time ordering. For a Gaussian pulse-antipulse sequence, one may show that 
V it = 2ahe~ 13 sin ^T s [a x sin 2jT — a y cos2'yT], where T — (T\ + T 2 )/2. Note that for 



t — T 2 ^> r, Vit depends on T\ and T 2 but not t. Then, e 



-iVitjn 



cos 



2ae * sin7T s /2 



a x sin 27T — a y cos 27T ■ sin 2ae @ sin ^T s . Consequently, 



u?(t) 



cos[2ae ^ sin7T s ] e 2i7T sin[2a;e ^ sin7T s 
— e 2ryT sin[2ae - ^ sin7T s ] cos[2ae _/3 sin7T s ] 



(26) 



We use this result in the next section to study the effect of time ordering on the transfer of 
population from one state to another. 

r i 



u°(t) 



We note that as either 7 



or a — > 0, one has Uj(t) 



1 



in contrast to 







-iyt 



. For simplicity, we may consider the case where each pulse is an 



ideal kick, i.e. (3 = 0. Then in the perturbative regime, expanding in a, one finds Uj n ~ 
I — 2a 2 sin 2 7T.,, which can be compared to e - * 7 *?/^ antlklck « 1 — 2ia 2 e~ nTs sin 7T S , so that 
except for special values of 7 (including 0), these matrix elements differ by 0(a 2 ). Similarly, 
U?io w 2asin^T s e- 2i ^ T + 0(a 3 ) compared to e'^U^ antikick w 2asin7T s e" 2 ^ + C(a 3 ). 



; /12 



Hence C// and Uj agree to leading order in a. This should not be surprising, as time ordering 
has no effect at leading order in perturbation theory. 
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III. CALCULATIONS 

As an illustrative specific example we present in this section the results of numerical 
calculations for 2s — ► 2p transitions in atomic hydrogen caused by a Gaussian pulse of width 
r. The occupation probabilities of the 2s and 2p states are evaluated by integrating two-state 
equations using a standard fourth order Runge-Kutta method. This enables us to verify the 
validity of our analytic solutions for kicked qubits in the limit r — ► and also to consider 
the effects of using finite-width pulses. In this system, the unperturbed level splitting is the 
Lamb shift, AE = E 2p — E 2s = 4.37 x 10 _6 eV. The corresponding time scale is the Rabi 
time, Tab = 2nh/AE = 972 x 10~ 12 s, which gives the period of oscillation between the 
states. 

For any practical system, the pulse duration r can neither be too large nor too small. If r 
is larger than T&e, then the pulse will not be sudden and the kicked approximation will fail. 
On the other hand if r is too small, then the interaction will have frequency components that 
couple the initial state to other levels. Specifically, if r is less than 2ith / (Ez p — E 2s ) ~ 1CT 15 s, 
then the interaction will induce transitions into the 2>p level and the system will not be well 
approximated by a two-state system. Also there is another constraint in our case. If the 
experiment lasts longer than the lifetime of the 2p state, 1.6 x 10™ 9 s, then we lose population 
from our two-state system, i.e. dissipation cannot be neglected. Similar calculations can be 



4]. 



done in many other applications, including, for example, Josephson junctions 

In the first part of this section we present results for the target state occupation proba- 
bility, P 2 , as a function of time. We shall examine how well the approximations we use are 
satisfied for a 2s — 2p transition caused by a pulse of finite width. We shall do this first for 
a single pulse and then for a double pulse. In the second part of this section we examine ef- 
fects of time ordering. Here we shall evaluate P 2 {t) both with and without time ordering for 
pulses of finite width. This will be done in both the Schrodinger and intermediate pictures. 

A. Pulsed two-state system 

In our numerical calculations we use for convenience an interaction of the form V(t) = 
(afi/y / 7rr)e _ ^ _Tft ) / T , i.e. a Gaussian pulse centered at T^ with width r. The evaluation 
of the integrated pulse strength a in terms of the dipole matrix element for the 2s — 2p 
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transition is discussed in a previous paper [18[. When r is small enough for the sudden, 
kicked approximation to hold, V(t) —>■ aM(£ — 2V), and the analytic expressions of Eqs. (JT7J1 
and (|23|) apply- Here we shall determine how the occupation probability P2W depends on 
the pulse width r, to find where the kicked results are approximately valid for finite pulses. 
We do this first for a single pulse chosen so that an ideal kick would transfer the occupation 
probability Pzif) suddenly from zero to one at t = T^. Then we consider two equal and 
opposite pulses occurring at times 7\ and T 2 . We study this doubly kicked system as a 
function of both pulse width r and separation interval T s — T 2 — T\. 

1. Single pulse 




100 150 200 
Time (ps) 



300 



FIG. 1: Occupation probability of the target state as a function of time for a qubit interacting with 
a single pulse. The heavy solid line corresponds to r = 1 ps (almost an ideal kick), the thin dashed 
line to r = 10 ps (where small deviations from an ideal kick occur), and the thin dotted line to 
r = 100 ps (where the kicked approximation is breaking down). The Rabi time for the oscillation 
between the states is 972 ps. 



In Fig. 1 we show results of a calculation for the probability Pi{t) that a hydrogen atom 
initially in the 2s state makes a transition into the 2p state when strongly perturbed by 
a single Gaussian pulse applied at t — T*.. We have obtained our results by numerically 
integrating the two-state coupled equations, 

1 



%hh\ 



a 



-AE ai + _ 

2 y/TTT 
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e -(t-n?/r\ 2 



iha 2 = -AEa 2 + _^ e -('- r *) 2 /- 2 ai 

2 y/lYT 



(27) 



Here the pulse is applied at Tk = 150 ps and we have chosen a = n/2 so that in the limit 
of a perfect kick all of the population will be transferred from the 2s to the 2p state after 
t = T k . 

In Fig. 1 one sees that the ideal kick results are very nearly achieved by choosing r to be 
a factor of 10 -3 times smaller than the Rabi time, Ta_e, in which the population oscillates 
between the 2s and 2p states. When t/T&e ~ 10~ 2 , a small deviation from an ideal kick can 
be seen in the figure. In this case PziTf) = 0.9977. When t/T^e ~ 10 _1 , the transition takes 
a few tenths of a nanosecond to occur and only 82% of the population is transferred at 300 
ps. The error in PziTf) resulting from the kicked approximation grows as (t/T^e) 2 ~ /3 2 , 
as expected from Eq. (|23|). 

2. A positive followed by a negative pulse 
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FIG. 2: Occupation probability of the target state as a function of time for a double pulse that 
returns the system to its initial state in the kicked limit. The heavy solid line corresponds to r = 1 
ps (almost ideal kicks), the thin dashed line to r = 10 ps (where small deviations from ideal kicks 
occur), and the thin dotted line to r = 100 ps (where the kicked approximation is breaking down). 

In Figs. 2 and 3 we show results of a calculation for the probability p2{t) that a hydrogen 
atom initially in the 2s state makes a transition into the 2p state when acted on by a double 
Gaussian pulse. Two Gaussian-form pulses are applied, the first at t — T\ and the second 
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FIG. 3: Occupation probability of the target state as a function of time for a double pulse that 
fully transfers population in the kicked limit. The heavy solid line corresponds to r = 1 ps (almost 
ideal kicks), the thin dashed line to r = 10 ps (where small deviations from ideal kicks occur), and 
the thin dotted line to r = 100 ps (where the kicked approximation is breaking down). 

at t — T 2 . The separation interval between pulses is T s = T 2 — T\. The final occupation 
probability of the target state is measured at t — Tf, The pulses are opposite in sign, but 
otherwise identical, so the interaction integrated over the whole interval [0, Tf] is zero. 
Our results for this double pulse have been obtained by numerically integrating 



iha\ = —-AEai + (a/y/nr) 
%ha2 = -AEa 2 + (a/y/nr) 



^) -r-c-^r 



't-Ti \2 



■t-T 2 - 



a 2 
Oi 



(28) 



In Fig. 2, the first pulse is applied at 7\ = 100 ps and the second at T 2 = 586 ps, giving a 
separation time T s = 486 ps = tt/2'j. From Eqs. (JHJ and (f23j) one sees that this is precisely 
the value of ^T s required to yield complete transfer from the 2s to the 2p state at T\ and then 
full transfer back to the 2s state at time T 2 in the limit of an ideal double kick. Moreover, 
we have chosen an integrated pulse strength a = tt/2 for each pulse for the same reason. 
The parameters for Fig. 3 are the same except that a = 7i~/4, so that the 2p target state is 
fully populated after T 2 for an ideally kicked system. 

As in Fig. 1 we see that when tJT^e < 10~ 3 , the kicked limit is well satisfied, and when 
t/T/^e ~ 10~ 2 , small deviations from an ideal kick can be seen in Fig. 2. In this case 
-^2(^7) — 11 x 10~ 5 . Again deviations from the ideal kick limit in the transfer probability 
P 2 grow as (r/T AE ) 2 . 



17 



The results in Fig. 3 are for a double pulse that first takes the population halfway from 2s 
to 2p, and then the rest of the way for an ideal kick-antikick sequence. When t/T&e ~ 10~ 2 , 
small deviations from an ideal kick can be observed in the figure near both steps. At 700 
ps, P 2 (Tf) = 0.99934, i.e. the population is nearly, but not quite perfectly transferred to 
the target state. When t/T^e ~ 10 -1 , the transition takes a few tenths of a nanosecond to 
occur and only 80% of the population is transferred, i.e. the transfer is not ideal. 

B. Time ordering 

In this subsection we consider the more complex issue of time ordering in 2s — 2p tran- 
sitions in atomic hydrogen caused by a single or double Gaussian pulse. Finite-width pulse 
effects are again considered. The effect of time ordering is evaluated by comparing results of 
calculations with and without time ordering for the probability P 2 of transferring an electron 
population from the 2s launch state to the 2p target state. Since the limit without time 
ordering is different in the Schrodinger and intermediate pictures, we include results for both 
pictures. 

The equations including time ordering are given by Eqs. ()27j) and (|2*B*j) above. The 
analogous equations without time ordering are found by taking V(t) — > V in the Schrodinger 
picture and Vi(t) -^Vi in the intermediate picture. 

1. Single pulse 

In Fig. 4 we show the effects of time ordering for a single pulse. We have chosen our 
parameters so that in the limit of an ideal kick the population is completely transferred 
from the 2s to the 2p state of hydrogen at time Tj,, as described above. On the left hand 
side of Fig. 4 we show how the occupation probability of the 2p target state varies as a 
function of the Gaussian pulse width r for three different values of the observation time Tf. 
For sharp pulses the exact transfer probability P2(Tf) and the transfer probability without 
time ordering in the intermediate picture Pj 2 (Tf) are quite similar, but differences appear, 
as expected, when t/T&e = /3/tt becomes large. 

However in the Schrodinger picture there are very large differences between the results 
with and without time ordering, P 2 (Tf) and P^Tf), even for an ideal kick. This occurs 
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FIG. 4: Target state probability as a function of the pulse width r (on the left), and as a function of 
the observation time Tf (on the right). Here Ta_e = 2Trh/AE = 71/7 is the Rabi time for oscillations 
between the states, where AE = E<i v — E<i s - The heavy line denotes probability including time 
ordering, the dashed line denotes the probability in the intermediate picture without time ordering, 
and the dotted line represents the probability in the Schrodinger picture without time ordering. 
On the right, the lines begin at the midpoint of the pulse, Tf = T&. The Schrodinger results damp 
out for large Tf on the right as explained in the text. 

because the energy splitting AE is non-zero, and for Tf > ah/ AE = oT^e/^-, the average 
potential V = ot/Tf becomes smaller than the energy splitting AE. Thus, for a given pulse, 
the influence of the potential necessarily decreases at large Tf, and any transfer probability 
becomes exponentially small. In effect, the free propagation before and after the pulse 
diminishes the effect of the pulse itself in the Schrodinger picture, when time ordering is 
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removed. This behavior contrasts with the intermediate picture result (Eq. (|llj)). where 
Pj 2 (Tf) depends on (3 = 7 r but not on Tf, as seen also on the left side of Fig. 4. The 
contrast is evident on the right hand side of Fig. 4 where, after the pulse has died off, the 
value of P 2 (Tf) dies out as Tf increases, while Pj 2 (Tf) approaches a constant. 

For a single narrow pulse one may compare the time-ordered result for the transfer prob- 
ability using the kicked approximation (Eqs. (|T7|) and (J22|)) with the exact expressions in 
the absence of time ordering in the Schrodinger and intermediate pictures, given by Eqs. (jHJ) 
and (JID), 

P 2 {T f ) = sm 2 a + 0(a 2 (3 2 ) 

PS(T f ) = a2 f m)2 sin 2 y^ + ( 7 T /)2 (29) 



PU?f) 



sin (ae 



These three equations are consistent with the numerical results shown in Fig. 4. As AE — > 0, 
jTf and (3 become small, time ordering effects disappear, and all three results coincide at 



sin 2 a. 



2. A positive followed by a negative pulse 

Finally, we consider the role of time ordering in the case of two equal and opposite pulses 
separated by an interval T s . In the limit of ideal kicks, the kick-antikick evolution operator 
has been expressed analytically above with and without time ordering in Eqs. (|24j) and (|26|). 
This yields for the probability P 2 (Tf) at times Tf after the second kick, 

P 2 (T f ) = sin 2 7 T s sin 2 2a + 0(/5 2 ) 
P? 2 (T f ) = sin 2 [2ae~P 2 sin 2 ( 7 T s )] (30) 

P°(T f ) = . 

The existence of these analytic results is helpful in studying the role of time ordering in P 2 . 
In Fig. 5 we compare P 2 with Pf 2 as a function of the separation time T s between the 
pulses. Here the measurement time Tf is taken to be well after the second pulse has decayed, 
Tf — T 2 ^> r. We note that, as expected, in all cases the occupation probability goes to zero 
as the two opposite pulses coalesce, i.e. as T s — > 0. On the left side we show the occupation 
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FIG. 5: Target state probability as a function of the separation time T s = T2 — T\ between the 
two pulses, for integrated pulse strength a = n/2, vr/4, and 3tt/8. On the left the pulse width is a 
rather narrow 10 ps, while on the right the pulse width is 100 ps, where the kicked approximation 
is breaking down. Again the heavy line denotes the exact result including time ordering, while the 
thin dashed line denotes the intermediate picture result without time ordering. The Schrodinger 
picture result without time ordering is identically zero for all values of T s and a. 

probability for a Gaussian pulse of width r = 10 ps for three values of the pulse strength a, 
namely a = tt/2 corresponding to a kick that turns a qubit from off to on at T\ and back 
to off at T 2 , a — tc/A where a qubit is on after T 2 , and an intermediate strength a = 37r/8. 
On the left side r /T& E = 10~ 2 , so that the kicked result is accurately obtained. We note 
that Pj 2 ~ P2, and time ordering effects vanish, when T s /T&e = ^T s /pi is an integer or 
half-integer, consistent with Eq. (|3U|) . Away from these special values, large time ordering 
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effects are present. On the right side, the pulses are quite broad (t/T^e = 10 x ), and the 
kicked approximation is clearly breaking down. 

IV. DISCUSSION 

In this paper we have studied the role of time ordering in a strongly perturbed two- 
state quantum system. We have defined the time ordering effect as the difference between 
a calculation with time ordering and one without time ordering. This is the way correla- 
tion, entanglement, and non-random processes are also defined J26(. In all of these cases 
it is useful to define carefully the limit without the effect (time ordering, correlation, or 
non-randomness). In the case of time ordering we have seen that the limit without time cor- 
relation depends on how one separates H into Hq + V, i.e. what we call a choice of gauge. In 
practice this choice sometimes rests on the choice of the time averaged interaction. Precisely 
this same issue occurs in defining correlation |27| . namely the somewhat arbitrar y ch oice 
of a mean field interaction. A similar problem can arise in defining entanglement [28] and 
non-random processes. In any case the problem is not new. 

For clarity and simplicity we have chosen Hq = const x o z and V = f(t)Vo<j x in this 
paper. However other gauge choices are possible and in some cases may be more sensible. 
Specifically, in many cases experimental conditions can lead to a sensible gauge choice where, 
for example, the asymptotic state of an unperturbed atom is an eigenstate of Hq, and V 
corresponds to an external electric or magnetic field imposed on the atom during part of the 
experiment. Such a V may sensibly contain a z and/or a y components. 

We have shown above that the limit of no time ordering depends on the picture (rep- 
resentation) used. While this difference is small when time ordering effects are small, the 
differences can be large otherwise. Hence it might be argued that our time ordering analysis 
is useful primarily to determine if time ordering effects are small or large. This argument 
sometimes occurs in the use of correlation (although it is not common to use different rep- 
resentations for correlation). In any case the dependence on representation is not new [29]. 
Convergence properties of the Magnus expansion in the Schrodinger and interaction pictures 
have been known for some time to differ widely [30]. While the Schrodinger picture result 
(that in a sense corresponds to an especially simple gauge choice of H Q = 0) is formally 
easier to write down, it often does not yield predictions as reliable as the more complete 



22 



results found in the intermediate picture. We have illustrated this above with calculations 
and analysis of U° and Uf for single and multiple pulses. The intermediate picture takes 
maximum advantage of knowledge about the eigenstates and spectrum of H . In other words 
the intermediate picture is generally more complete than the Schrodinger picture and often 
more sensible. It is useful to separate H into H + V in such a way as to include as much 
of the problem as possible in Ho, whose solutions are known. In the extreme limit of the 
Heisenberg picture (which can be thought of as a gauge choice where H = Hq and V = 0), 
we have Uh = U^ = 1, and there is never any time ordering in the time evolution. 

Here we have worked in the time domain and formulated the question of time ordering by 
explicitly working with the Dyson time ordering operator, T. Equivalently one may work in 
the energy domain, as has been done recently in the context of atomic scattering to analyze 
experimental data and identify time ordering effects 
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14 . A key transformation for time 



ordering from the time domain to the energy domain is the Fourier transform of the step 
function, namely, / e tE,t O(t)e~ lEt dt = E _^, +i = 7t5(E — E') + iP v ^^j, where the effect 
of time ordering is associated with irj, which gives rise to the principal value term. The irj 
carries the effect of the boundary condition on the Green's function in energy space. This is 
discussed in more detail elsewhere |8lll0|. Since \H{t"), H{t')] provides a connection between 
interactions at t" and t', the quantum time propagator includes non-local effects in time. 
An example of a counter-intuitive time sequence occurs in stimulated Raman adiabatic 
passage (STIRAP), where efficient and robu st p opulation transfer is attained using two 
pulsed radiation fields in a three-level system 

The time ordering effects considered in this paper are associated with a sequential ordering 
of interactions. The normal boundary condition imposed on the evolution operator is that 
the sequence proceeds in the direction of increasing time (or alternatively decreasing time 
to study time reversal). Hence effects of time ordering are associated with a direction of the 
flow of time. In this paper we have not included dissipation; hence all amplitudes explicitly 
satisfy invariance under time reversal, e.g. in Eqs. (fT7j) and (J23J) which include effects due 
to time ordering. Effects of time ordering have been observed in systems that satisfy time 
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reversal invariance 
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This means that observable evidence of the direction of time 



(time ordering) can be obtained without violating the symmetry of time reversal invariance. 
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V. SUMMARY 

In this paper we have given a definition of time ordering in a strongly perturbed quantum 
system, namely that time ordering is the difference between calculations with and without 
time ordering. This definition is similar to the definition of correlation. In both cases effects 
arise from differences between an instantaneous interaction and its averaged value. When the 
effect of time ordering is small, the dependence on representation is weak. However, when 
time ordering effects are large, the difference between representations can also be large. We 
have considered in detail time ordering for qubits that are strongly and suddenly perturbed 
by an external interaction. We have illustrated our methods for a 2s — 2p transition in atomic 
hydrogen caused by a Gaussian pulse of finite width in time. Other diabatically changing 
qubits may also be analyzed with our methods. Simple analytic expressions have been 
given for the occupation amplitudes and probabilities for kicked qubits, including single and 
multiple kicks. We think that it should be possible to find analytic solutions for correlated 
kicked qubits, so that time coupled interacting qubits may also be studied analytically. 
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